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CN ■ Abstract 

We propose a generalisation of the Faddeev-Popov trick for Yang-Mills fields in 
the Landau gauge. The gauge-fixing is achieved as a genuine change of variables. In 
particular the Jacobian that appears is the modulus of the standard Faddeev-Popov 
determinant. We give a path integral representation of this in terms of auxiliary 
bosonic and Grassman fields extended beyond the usual set for standard Landau 
gauge BRST. The gauge-fixing Lagrangian density appearing in this context is local 
and enjoys a new extended BRST and anti-BRST symmetry though the gauge-fixing 
0^ ■ Lagrangian density in this case is not BRST exact. 
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^ ■ 1 Introduction 



The elevation of Faddeev-Popov (FP) gauge-fixing of Yang-Mills theory be- 
yond the realm of perturbation theory has been intensely pursued in recent 
years for many reasons. Nonperturbative gauge-fixed calculations on the lat- 
tice are being compared to analogous solutions of Schwinger-Dyson equations 
[1,2]. As well, the long-term goal of simulating the full Standard Model using 
lattice Monte Carlo requires the Ward-Takahashi identities associated with 
BRST symmetry [3] in order to control the lattice renormalisation. The main 
impediment to nonperturbative gauge-fixing is the famous Gribov ambiguity 
[4]: gauges such as Landau and Coulomb gauge do not yield unique represen- 
tatives on gauge-orbits once large scale field fluctuations are permitted. To 
some extent one could live with such non-uniqueness if one could incorporate 
all Gribov copies in a computation. However the no-go theorem of Neuberger 
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[5] obstructs even this: (a naive generalisation of) BRST symmetry forces a 
complete cancellation of all Gribov copies in BRST invariant observables giv- 
ing 0/0 for expectation values. In particular, Gribov regions contribute with 
alternating sign of the FP determinant. 

Here we shall propose an approach which takes seriously that gauge-fixing 
when seen as a change of variables involves a Jacobian being the absolute value 
of the Faddeev- Popov determinant. Usually the absolute value is dropped 
either because of an a priori restriction to perturbation theory or because of 
the identification of the determinant in terms of an invariant of a topological 
quantum field theory [6] such as the Euler character [7,8]. In the latter case 
the Neuberger problem is encountered. 

The approach we describe in the following is not restricted to perturbation 
theory. Moreover, because it will be seen to involve a gauge-fixing Lagrangian 
density that is not BRST exact it falls outside the scope of the preconditions 
for the Neuberger problem. In the next section we shall derive the Jacobian 
associated with gauge-fixing in the presence of Gribov copies. We shall give 
a representation of the "insertion of the identity" in this case in terms of a 
functional integral over an enlarged set of scalar and ghost fields. The extended 
BRST symmetry of this new gauge-fixing Lagrangian density will be described 
though we will see that the final form of the gauge-fixing Lagrangian is not 
BRST exact. 



2 Field theoretic representation for the Jacobian of FP gauge fixing 



In the following we shall formulate the problem in the continuum approach to 
gauge theory. 

Our aim is to generalise the standard formula from calculus for a change of 
variable: 

-i 

= f d Xl ...dxJ^(f(x)). (1) 

f=0 




Here one is changing from integration variables x to those satisfying the con- 
dition f(x) = and where, for Eq.(l) to be valid, in the domain of integration 
of x there must be only one such solution. In the context of gauge-fixing of 
Yang-Mills theory the generalisation of Eq. (1) is 

-i 

= fvgS[F[«A]] (2) 

F=0 
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where represents the gauge field, g is an element of the SU(N) gauge 
group, Vg is the functional integration measure in the group and 



F[ 9 A] = 



(3) 



is the gauge- fixing condition. We shall be interested in Landau gauge F[A] = 
d^A^. As in the calculus formula, here Eq.(2) is only valid as long as Eq.(3) 
has a unique solution. This is known not to be the case for Landau gauge. The 
FP operator nevertheless is MpL4] = (5F[ 9 A]/5g)\F=o and its determinant is 
A F [A] = det(M F ). For the Landau gauge M F [A] ab = d^Df[A\ with Df[A\ 
the covariant derivative with respect to A® in the adjoint representation. Now 
the standard FP trick is the insertion of unity in the measure of the generating 
functional of Yang-Mills theory realised via the identity (which follows from 
the above definitions): 



By analogy with standard calculus, in the presence of multiple solutions to 
the gauge-fixing condition Eq.(4) must be replaced by 



where iV F L4] is the number of different solutions for the gauge-fixing condition 
F[ 9 A] = on the orbit characterised by A, where A is any configuration on 
the gauge orbit in question for which det M F ^ 0. 

It is known that Landau gauge has a fundamental modular region (FMR), 
namely a set of unique representatives of every gauge orbit which is moreover 
convex and bounded in every direction [9,10]. The following discussion can 
be found in more detail in [11]. Denoted A, the FMR is defined as the set of 
absolute minima of the functional Va[<?] = / d^x^A) 2 with respect to gauge 
transformations g. The stationary points of Va[<?] are those A^ satisfying the 
Landau gauge condition. The boundary of the FMR, <9A, is the set of degen- 
erate absolute minima of Va[<7]. A lies within the Gribov region Qq where the 
FP operator is positive definite. The Gribov region is comprised of all of the 
local minima of V^[g]. The boundary of Q , the Gribov horizon dflo, is where 
the FP operator M F (which corresponds to the second order variation of Va[<?] 
with respect to infinitesimal g) acquires zero modes. When the degenerate ab- 
solute minima of dA coalesce, flat directions develop and M F develops zero 
modes. Such orbits cross the intersection of d A and 8Qq. The interior of the 
fundamental modular region is a smooth differentiable and everywhere con- 
vex manifold. Orbits crossing the boundary of the FMR on the other hand 




(4) 




(5) 
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will cross that boundary again at least once corresponding to the degenerate 
absolute minima. 



Though, at present, there is no practical computational algorithm for con- 
structing the FMR, it exists and we will make use of it for labelling orbits, 
i.e., A u are defined to be configurations in the FMR, A u e A. Since every 
orbit crosses the fundamental modular region once we are guaranteed to have 
N F > 1. In turn the 9 A U fulfilling the constraint of Eq. (3) would be every 
other gauge copy of A u along its orbit. Eq.(5) is equal to the number of Gribov 
copies on a given orbit, Nqc — N F — 1, except that copies lying on any of the 
Gribov horizons (A F = 0) do not contribute to N F . 

The finiteness of N F in the presence of a regularisation leading to a finite 
number of degrees of freedom (such as a lattice formulation) can be argued 
as follows. Consider two neighboring Gribov copies corresponding to a single 
orbit. If they contribute to N F they cannot lie on the Gribov horizon. Therefore 
they do not lie infinitesimally close to each other along a flat direction, namely 
they have a finite separation. This is true then for all copies on an orbit 
contributing to N F : all copies contributing to N F have a finite separation. 
But the g which create the copies of A n belong to SU(N) which has a finite 
group volume. Thus for each space-time point there is a finite number of such 
g. We conclude then for a regularised formulation that N F is finite. 

Consider then the computation of the expectation value of a gauge-invariant 
operator 0[A\ over an ensemble of gauge- field configurations A u which is this 
set of unique representatives of gauge orbits discussed above. 

Note that for a gauge-invariant observable, it makes no difference whether 
A u G A or if the A u J s are any other unique representatives of the orbits. 

The expectation value on these configurations 



J VA u O[A u ]e-^ 
JVA u e- 



(o[a\> = j z:z YM ( 6 ) 



is well-defined. Since in any regularised formulation N F is a finite positive 
integer, we can legitimately use Eq.(5) to resolve the identity analogous to the 
FP trick and insert into the measure of integration for an operator expectation 
value. We thus have 

{Q[A]) = SVAv^SVgSjFpA]) | det M F \°A]\ Q[A] e~s™W ^ 



JVA ujr ^ x - l JVg5(FM) det M F [sA] e s™lA] 

We can now pass A?^[A U ] under the group integration Vg and combine the 
latter with VA U to obtain the full measure of all gauge fields V( 9 A U ) which 
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we can write now as DA. N F is certainly gauge-invariant: it is a property of 
the orbit itself. So N F [A U ] = N F [ 9 A U ] = N F [A\. Thus we can write 



(0[A\) = 





det M F [°A] 


0\A] e -^M[A] 


fVAj^5(F[9A]) 


det M F [aA] 


q—Sym[A] 



(8) 



Perturbation theory can be recovered from this of course by observing that 
only A fields near the trivial orbit, containing A = and for which 5Vm[/4] = 0, 
contribute significantly in the perturbative regime: the curvature of the orbits 
in this region is small so that the different orbits in the vicinity of A = 
intersect the gauge-fixing hypersurface F = the same number of times. Then 
the number of Gribov copies is the same for each orbit, N F is independent of 
A u and we can cancel N F out of the expectation value. In that case 



(0[A\) = 



JVA5(F[A}) 


det M F [A\ 


0[A] e - SYM ^ 


JVA5(F[A}) 


det M F [A] 





(9) 



In turn, observing that fluctuations near the trivial orbit cannot change the 
sign of the determinant, the modulus can also be dropped and one recov- 
ers the usual starting point for a standard BRST invariant formulation of 
Landau gauge perturbation theory. Note that perturbation theory is built on 
the gauge-fixing surface in the neighbourhood of A = 0, which for a gauge- 
invariant quantity will be equivalent to averaging over the Gribov copies of 
A = as in Eq. (9). For the non-perturbative regime, the orbit curvature in- 
creases significantly and in general there is no reason to expect that N F would 
be the same for each orbit. Moreover the determinant can change sign. 

Let us focus on the partition function appearing in Eq. (8) 

^gauge-fixed = / VA N F 1 [A] \ det (M F [A] ) 1 6(F[A] ) e^" (10) 



The objective is to generalise the BRST formulation of Eq.(10) such that it 
is valid beyond perturbation theory taking into account the modulus of the 
determinant. We thus start with the following representation: 

det(M F L4]) = sgn(det(M F L4])) det(M F L4]). (11) 



As mentioned, the factor det(M F [A|) in Eq.(ll) is represented as a functional 
integral via the usual Lie algebra valued ghost and anti-ghost fields in the 
adjoint representation of SU(N). Let us label these as c a ,c a . It is usual also 
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(see for example [12]) to introduce a Nakanishi-Lautrup auxiliary field b a . Thus 
the effective gauge-fixing Lagrangian density 

£dct = -Wd^Al + h) a b a + c a M F b c b (12) 



yields [12] 

limy Vc a Vc a Vb a e-J d4xCd * = S(F[A}) det(M F [A}). (13) 



In order to write the factor sgn(det(MirL4])) in terms of a functional integral 
weighted by a local action, we consider the following Lagrangian density 

£ sgn = iB a M ab ^ h - id a M ab d b + l -B a B h (14) 



with d a , d a being new Lie algebra valued Grassmann fields and f a , B a being 
new auxiliary commuting fields. Consider in Euclidean space the path integral 

Z sgn = J Vd a Vd a V(p a VB a e-f d4xCs * n . (15) 



Completing the square in the Lagrangian density of Eq.(14), the B field can 
be integrated out in the partition function leaving an effective Lagrangian 
density 

4gn = ^((M f ) r ) flt MV - id a Mfd b , (16) 



where (Mp) T denotes the transpose of the FP operator. Integrating all remain- 
ing fields now it is straightforward to see that the partition function Eq.(15) 
amounts to just 



Z sgQ = d6t(MF) = sgn(det(M F )) . (17) 
y/det((M F ) T M F ) 

Thus the representation Eq.(15) can be used for the first factor of Eq.(ll). The 
Lagrangian density of Eq.(14) therefore combines with the standard BRST 
structures of Eq.(12) coming from the determinant itself in Eq.(ll) so that an 
equivalent representation for the partition function based on Eq.(lO) is 



^gaugc 



-fixed = f VA a ^c a Vc a Vd a Vd a Vb a Vcp a (N F [A\)- 1 e- SYM - Sd <' t - Sa ^(18) 
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with Sdet and S sgn the actions corresponding to the above Lagrangian densities 
Eqs. (12,14). 



3 A new extended BRST 



The symmetries of the new Lagrangian density, C sgn , are essentially a boson- 
fermion supersymmetry and can be seen from Eq.(14). In analogy to the stan- 
dard BRST transformations typically denoted by s, we shall denote them by 
the Grassmann graded operator t 



tip a = d a 
td a = 
td a = B a 
tB a = 0, 



such that 



tC =0 

l " < — sgn u 



(19) 



(20) 



and trivially t£ Y M = 0. Eqs. (19) realise the infinitesimal form of shifts in the 
fields. The operation t is nilpotent: t 2 = 0. Using Eqs. (19) we can give the 
following form for the Lagrangian density £ sgn , 



r 

sgn 



t (d a {iMfy h + ^B a )^j . (21) 



The question now is how to combine this with the standard BRST transfor- 
mations 



sa; = D; b c b 

sc a = -^gf abc c b c c 
sc a = b a 

sb a = 0. (22) 



The transformations due to t and s are completely decoupled except that the 
latter also act on the gauge field on which the FP operator Mp depends. We 
propose the following unification of these symmetry operations. Consider an 
operation S block- diagonal in s and t: S — diag(s, t). The operator acts on the 
following multiplet fields: 
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A a 



( A a 



c a 



d a 



, C° 



( c a 



V 



d a 



, B a 



b a 



(23) 



We see that these fields transform under S completely analogously to the 
standard BRST operations 



SA a = V ab C b 

Qna ^nabcnbnc 

SC a = B a 

SB a = 0, (24) 
where k — 1, 2 label the elements of the multiplets, and 



V ab = dieLg{Df ,5 ab ' 
1 
2 



(25) 



Note that nilpotency is satisfied, S 2 = 0. We shall refer to this type of op- 
eration as an extended BRST transformation which we distinguish from the 
BRST-anti-BRST or double BRST algebra of the Curci-Ferrari model [13,14]. 
We can thus formulate the gauge-fixing Lagrangian density for the Landau 
gauge as 



Tr S 



^ c a F a 







\ 



v d a {iMf^ b + \B a ) / 



(26) 



This approach admits also an extended anti-BRST operation: 



SA a = V ab C b 

Qna ^abcnbnc 

SC a = -B a 
SB a = 0. 



(27) 



Writing S = diag(s, t) we can extract the standard anti-BRST s-operations 
[14,15] 



SAI = Dfc b 
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sc a = - l -gt b ^c c 
sc a = -b a 

sb a = (28) 
and those corresponding to t: 

iip a = d a 
id a = 
id a = -B a 

iB a = 0. (29) 



Moreover, the ghosts and anti-ghosts in this extended structure also fulfill the 
criteria for being Maurer-Cartan one-forms, 



SC + SC = 0. 



(30) 



However there is no extended BRST-anti-BRST (or double) symmetric form 
of the gauge- fixing Lagrangian density Eq. (26), unlike the two pieces of which 
it consists. Such a representation exists in the s— sector of Landau gauge: 



C gi , s = -ssA^. 



(31) 



In the t— sector, the corresponding structure is 

1 _ 



£gf,t — -tt 



tp a M£ b V b + d a d a 



(32) 



However the complete Landau gauge-fixing Lagrangian density can only be 
expressed via a trace, namely as 

£ gf = ^TrSSW (33) 



with 

W = diag (A^Al, ip a M£ b ip b + d a d a ) . (34) 

Nevertheless this compact representation formulates the modulus of the de- 
terminant in Landau gauge fixing in terms of a local Lagrangian density and 
follows as closely as possible the standard BRST formulation without the 
modulus. 
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4 Discussion and Conclusions 



We have thus found a representation for Landau gauge-fixing corresponding to 
the FP trick being an actual change of variables with appropriate determinant. 
The resulting gauge-fixing Lagrangian density enjoys a larger extended BRST 
and anti-BRST symmetry. However it cannot be represented rigorously as 
a BRST exact object, rather the sum of two such objects corresponding to 
different BRST operations. This means that some of the BRST machinery is 
not available to this formulation, such as the Kugo-Ojima criterion for selecting 
physical states. We discuss cursorily now the perturbative renormalisability of 
the present formulation of the theory. Note that the procedure leading to Eq. 
(26) does not introduce any new coupling constants; only the strong coupling 
constant g is present in MpL4] coupling the Yang-Mills field to both the new 
ghosts and scalars. The dimensions of the new fields are 

[ V ]=L\ [d] = [d]=L-\ [B]=L- 2 . (35) 

Most importantly in this context, the kinetic term for the new boson fields ip a 
is quartic in derivatives: 

A™ = y? a (d 2 )V, (36) 

which is renormalisable, by power counting, since ip a are dimensionless. Such a 
contribution is seemingly harmless in the ultraviolet regime: for large momenta 
propagators will vanish like 1/p 4 . Moreover it should play an important role in 
guaranteeing the decoupling of such contributions in perturbative diagrams. 
That such a decoupling should occur is clear from Eq. (11): in the perturbative 
regime fluctuations about = will not feel the sgn(detMp[A]), so that the 
field theory constructed in this way must be equivalent to the perturbatively 
renormalisable Landau gauge fixed theory. For example in the computation 
of the running coupling constant we expect that this property will lead to a 
complete decoupling of the t-degrees of freedom so that the known Landau 
gauge result emerges from just the gluon and standard ghost sectors. Naturally, 
the new degrees of freedom will be relevant in the infra-red regime, which will 
be the object of future study. 
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